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Abstract. Modeling important engineering problems related to flow-induced damage (in the context of 
hydraulic fracturing among others) depends critically on characterizing the interaction of porous media 
and interstitial fluid flow. This work presents a new formulation for incorporating the effects of pore 
pressure in a nonlocal representation of solid mechanics. The result is a framework for modeling fluid- 
structure interaction problems with the discontinuity capturing advantages of an integral based formulation. 
A number of numerical examples are used to show that the proposed formulation can be applied to measure 
the effect of leak-off during hydraulic fracturing as well as modeling consolidation of fluid saturated rock 
and surface subsidence caused by fluid extraction from a geologic reservoir. The formulation incorporates 
the effect of pore pressure in the constitutive description of the porous material in a way that is appropriate 
for nonlinear materials, easily implemented in existing codes, straightforward in its evaluation (no history 
dependence), and justifiable from first principles. A mixture theory approach is used (deviating only slightly 
where necessary) to motivate an alteration to the peridynamic pressure term based on the fluid pore pressure. 
The resulting formulation has a number of similarities to the effective stress principle developed by Terzaghi 
and Biot and close correspondence is shown between the proposed method and the classical effective stress 
principle. 



1. INTRODUCTION 

The use of hydraulic fracturing to increase yields from natural gas reservoirs has gained popularity re- 
cently due to changes in the economic climate, the desire to decrease dependence on foreign energy sources, 
and the fact that natural gas burns cleaner than coal when used to produce electricity. This renewed in- 
terest in hydraulic fracturing has presented a number of challenges associated with modeling the processes 
involved. In particular, the need exists to better model the evolution of damage in the reservoir to prevent 
unintentional penetration of underground water resources. This type of modeling also helps to establish 
injection parameters that will produce the best results and is useful in estimating yields from a particular 
well. Towards these modeling needs, this work presents a modified formulation for state-based peridynam- 
ics that includes the influence of fluid pore pressure in a form similar to the familiar principle of effective 
stress. Using the theory of interacting continua [TJ|2J[31II] we show that the principle of correspondence (the 
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idea that the state-based peridynamic representation of the constitutive law is equivalent to the classical 
rendition) holds for the proposed formulation. This is an important result in that the proposed formulation 
preserves the material behavior of the porous media as represented in the classical theory and provides a 
way to experimentally obtain the material parameters for the proposed formulation in a manner consistent 
with the classical approach. 

The development of the peridynamic theory in general is presented in [5, 6, [8] and has been applied 
to a number of engineering problems related to membranes, reinforced concrete, and carbon fiber nanotube 
reinforced composites. This work builds on the state-based (as opposed to bond-based) version of the theory 
that is presented in [BJ. The major advantages of the state-based approach include a material response 
that depends on collective quantities like volume change or shear angle and the ability to incorporate a 
constitutive model from the conventional theory of solid mechanics. The core contribution of this work is 
a modification of the force scalar state which introduces the effects of fluid pore pressure. We show, under 
certain conditions, that the proposed formulation is equivalent to the notion of effective stress as presented 
by Terzaghi and Biot O QUI [11], but when generalized can be used to incorporate other models for the 
interaction of the fluid and the porous media. In the present work we use an interaction model that is based 
on the drag produced by the fluid on the porous matrix, although the use of other models is equally feasible. 
To simplify the presentation, we assume that the fluid pressure has previously been evaluated numerically 
or analytically through some other means. In this way we avoid delving into the numerous complexities 
regarding a fully coupled analysis in which the fluid pressure is treated as an unknown. Although we make 
no mention of the stability of various coupling algorithms for this formulation, we intend to pursue this in a 
forthcoming work. 

A general overview of the theory behind poroelasticity is presented in [T5]. In subsurface geomechanical 
modeling poroelasticity, or the interaction of the pore fluid with the porous media plays an important role 
and has been the focus of several studies O HH E3 HBJ HZ1 HE] • In the present work we demonstrate the 
proposed method for a number of numerical examples motivated by problems found in [THl [301 UH EHj ■ In 
many cases analytic solutions are known for the particular problem. For the examples that do not have an 
easily derived analytic solution, the results of the proposed method are compared with a commercial finite 
element code in order to gauge accuracy. 

An outline of the proceeding sections is as follows: We begin with the development of the modified state- 
based peridynamic theory and its numerical implementation. Then, we present the governing equations 
in their discrete form. To demonstrate the interesting features of the proposed method we then show its 
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performance for several numerical examples. We then conclude with some general comments and discussion 
of the results. 



2. FORMULATION 

It will be helpful to begin with some basic definitions. The definition of a scalar state and a vector state 
are best understood in terms of their respective action upon a vector, The scalar valued image of the 
vector, £ under scalar state s is written: s(£). Likewise the vector valued image of the vector £ under the 
vector sate V is written V_(£). In general a state of order m is a function A(-) : T~L — > C m where H is 
a spherical neighborhood centered at the origin in M 3 and C m is the set of all tensors of order m. The 
governing solid mechanics equations, as posed in the state-based peridynamic theory, make frequent use of 
the definitions above. In this work we restrict our attention to quasi-static analysis and neglect the inertial 
term (although the formulation does not preclude its inclusion) . Consider a neighborhood, % x , of a particle 
with position, x (see Figure [T]). The state based equation of motion can be written in terms of a force vector 
state field, T that defines the collective force relationship of the particle at x with each neighbor at position 
x'. 



(1) 



/ {T[x,t](x' -x) -T[x',t](x-x')}dV x > +b(x,t) 







The radius of the spherical neighborhood H x , referred to as the horizon, will be denoted <5. The geometry 
of deformation is shown in Figure [2j Regarding the integration in the equation above, in this formulation 
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Figure 1. The spherical neighborhood, H x , and horizon, 6, of a particle at position x (the 
white particle at the center of the figure). The force interaction of this particle with the 
gray particles is considered non-negligible. 



we assume that all functions are Riemann integrable. The relative position of two particles in the reference 
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Figure 2. The geometry of deformation for peridynamics, (left) reference and deformed 
particle states (right) ordinary state based material response. 



configuration is denoted as £ 



(2) 



£ = x' — x 



Let y(x,t) represent the the deformed position of the particle originally at x 



(3) 



y(x, t) = x + u(x, t) 



where u(x,t) is the displacement. The reference position and deformation vector state fields, X and Y_ 
respectively, are defined as 



(4) 
(5) 



Y[x,t]{$ =y(x + t,t)-y(x,t) 



2.1. Linear elastic peridynamic solid. For an ordinary material, there exists a scalar state, t, acting in 
the direction of Y_ such that the following relationship holds. 



(G) 



T = t M 



where M(Y)(x' — x) is a unit vector pointing from the deformed position of x to the deformed position of 
x' . The force scalar state for a linear elastic peridynamic material is given by 



(7) 



t = — -lo x + Pw e d 
m 



where p = —kd is the peridynamic pressure, x is the reference position scalar state field, e d is the deviatoric 
extension scalar state field, u is an influence function (which for this work we assume is of unit value), m is 
the weighted volume, 9 is the dilatation, and k and j3 = 15 fi/m are material parameters. Via the principle 
of correspondence, it can be shown that k = represents the bulk modulus from the classical theory and 
fi = [A^ represents the shear modulus. 

The reference position scalar state field, x, is simply the magnitude of the reference position vector state 
field. Likewise, let y represent the magnitude of Y\ The deviatoric extension scalar state field, e d , stems 
from the extension scalar state field, e = y — x, as follows 

(8) , d = ,- 9 -§ 

2.2. Effective stress principle. Thus far we have made no deviation from the formulation presented 
in [B]) for an ordinary linear elastic peridynamic material. We now present a new peridynamic material, 
representative of a saturated porous medium effectively loaded by the fluid pressure in the pores. Since the 
force scalar state in the standard linear elastic peridynamic material above has been split into a dilatational 
or volume changing term and a deviatoric term, the extension to a peridynamic material with effective stress 
is straightforward. The influence of the fluid pore pressure can be included by modifying the definition of 
the peridynamic pressure as follows 

(9) p = —k9 + jpf 

where pf represents the fluid pressure and 7 is the fluid pressure coefficient which will be defined in a 
subsequent section. Note that the sign convention for the peridynamic pressure is positive in tension, while 
the convention for the fluid pressure is positive when it induces expansion of the bulk material. 

3. PERIDYNAMIC MATERIALS WITH EFFECTIVE STRESS 

In the present modification of the peridynamic material above, the fluid itself is not represented directly 
in the model, only its effect on the constitutive behavior of the porous medium. For this reason we assume 
that the fluid pressure is known a priori (through alternative means of analysis for the flow quantities). It 
is important to recognize that the effective stress principle does not alter the constitutive behavior of the 
material influenced by the pore pressure. Even though the loads produced by the fluid pressure are introduced 
in the constitutive relationship, they can easily be extracted to appear in the formulation as an applied 
load. We chose to introduce the pore pressure in the constitutive relationship to simplify the numerical 

implementation since no method for application of natural boundary conditions has been developed yet for 
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the peridynaimc theory. The saturated material under consideration behaves identically to its dry correlary 
as parameterized by its material constants, k and fj,. For this reason, the principle of correspondence holds 
(the idea that a peridynamic material is an equivalent representation of a material in the classical theory). 
To this effect, the formulation above does not preclude application to elastic-plastic materials, even though 
it has been presented here for a linear material to simplify the formulation. 

The discussion above suggests that to preserve correspondence and the constitutive behavior of the peri- 
dynamic material, 7 must be defined in such a way that its origins stem from the forces induced by the 
interaction of the fluid with the porous material. To derive 7 in such fashion we turn to the theory of 
interacting continua or mixture theory. 



3.1. Theory of interacting continua. The theory of interacting continua provides a method for char- 
acterizing the bulk response of a volume by incorporating the influence of the constituents of which it is 
composed. In the present formulation we use the theory to define an interaction force between the pore 
fluid and the solid matrix such that coupling takes place at the interface between them, but the effect is 
represented as a volumetric quantity. To present the theory of interacting continua in a consistent manner, 
we begin with a classical representation of the pore fluid and the porous media such that the porous media 
is treated as the classical equivalent of the ordinary linear elastic peridynamic material. At a later stage, we 
return to the peridynamic representation. 

In the theory of mixtures, the total density of the mixture and the mixture velocity are defined as: 



(10) 




where N is the total number of constituents. We assume the existence of a partial traction vector, *W and 
partial stress tensor, T^ l \ such that tW = T^n, where n is the normal vector to the surface. The total 



traction and total stress are written 



N 



N 



(11) 



t 




i=l 



T = 




The balance of mass for each constituent is enforced via 



(12) 



£)(i)p(i) 



+ p«div v (i) =m (i) 



Dt 
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where D D ^ is the total time derivative operator and the mass production terms from each constituent must 
balance (i.e. y\ —1 = 0). The balance of linear momentum for the individual constituents is written 



(13) 



div [t w 1 +c« =0 



In the above expression, c is the interaction term that represents the influence of the other constituents. 
Consider the case of two constituents, fluid flowing through a porous solid. We shall employ an interaction 
model of the following form 



(14) 



c = a(> (/) -t)W) 



The interaction term can be thought of as a drag term relative to the difference in the velocity of the fluid and 
the solid, where the interaction coefficient is a. The model for how the fluid and solid interact can be much 
more complex, for example we could consider localized changes in porosity resulting form the pore pressure, 
but for this work we consider only the simple case of a drag-like term. To enforce a Newtonian balance of 
forces the interaction terms for the fluid and the solid will be of opposite sign (i.e. = — = c). 

For the case of a simple fluid without distortional stresses, the stress in the fluid can be written T^' = 
—p^I. Note that I is the identity tensor (no italics). If we assume that the velocity of the solid is much 
smaller than the velocity of the fluid, t)W can be neglected in the interaction term and the balance of linear 
momentum for the fluid becomes Darcy's law [22] . 



(15) 



av U) _|_ g r ad p 



(/)! = 



The balance of mass for the fluid can be written 

d<t>^ 



(16) 



dt 



div 



,(/) 



77! 



(/) 



For a linear elastic solid, in the classical theory the stress tensor takes the form 

T 



(17) 



T (s) = A (s) div m (s) I + ^ (s) grad 



u 



grad 



u 



(<0 



where and fjf- 8 ^ are the material parameters (recall that the relationship between the bulk modulus 
and A^ s ) and fj,^ is = \^ + 2// s )/3). The balance of linear momentum for the solid is written 



(18) 



-air 



(«~div 
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In the absence of a fluid body force, equation ( 15 1 can be solved for the fluid velocity, av(f' = — grad \pW] ■ 
This result can be substituted into equation ( 18 ) to produce a modified form of the momentum balance 
expression for the solid 



(19) 



+ p (s) 6 (s) =0 



where "* = — p^I. In the expression above we have made use of the identity grad[p^] = div[p(^I]. 

In the peridynamic theory, direct application of the above analysis implies that the fluid pressure coeffi- 
cient, 7, in equation ^ is 1.0, which represents exactly Terzaghi's notion of effective stress. It is well known 
that a unit value for the pressure coefficient in Terzaghi's method performs well for granular materials with 
very small ratios of the parameter, K, for the bulk material versus the solid constituent alone. In order to 
apply the proposed formulation effectively to materials like stone or marble, we introduce a pressure coef- 
ficient that represents Biot's concept of effective stress, whereby the ratio of the moduli between the bulk 
material and the solid constituent alone is taken into account 

(20) 7 = 1- h 



where K is the bulk modulus of the volumetric material and K^ s > is the bulk modulus of the solid constituent. 

Remark 1. It is important to recognize that material behavior itself is unmodified. We merely chose to in- 
troduce the effect of the fluid pore pressure by means of a modified force scalar state, therefore correspondence 
still holds for our original ordinary linear elastic peridynamic material. 

3.2. Discrete form of the governing equations. The discrete form of the governing integral equation is 
written 

N 

(21) Yl *] « - x )- 2>i> - <)} AV < +b(x,t)=0 

i=0 

where N is the number of neighbors in the spherical neighborhood of x and AV X >. is the volume of cell 
representing the particle at x\ . The vector force state T, is evaluated in terms of the accumulated deformation 
of all the particles in the neighborhood of x according to equation ^ . 

4. NUMERICAL RESULTS 

This section presents a number of numerical examples that demonstrate the utility of the proposed method. 
To simplify the presentation, for all of the examples a particular fluid pressure field is prescribed as if the 
flow characteristics for the given problem are known. In a forthcoming work the coupling of the proposed 



method with existing flow modeling frameworks is investigated in detail. Since the formulation is inherently 
three-dimensional, all of the one-dimensional verification problems were performed on a domain of unit 
thickness. The exact solutions to the verification problems were derived using partial differential equations. 
Although the proposed method is integral based, the results should match under certain conditions due to 
the principle of correspondence. 



4.1. Lighthouse consolidation. The proposed method can be applied to consolidation of fluid saturated 
geologic material. The lighthouse consolidation problem represents a one-dimensional steady-state fluid- 
structure interaction problem with a known analytic solution. This problem is motivated by an example 
presented in |21j , and entails determining the deformation response of a partially submerged column of rock 
supporting a lighthouse. Including the influence of effective stress, the portion of the rock column below the 
water level experiences less deformation due to the pressure of the fluid in the pores. Above the water level, 
the rock column experiences increased subsidence due to capillary effects. For a column with cross-sectional 
area, A, length a above the water level, length b below the water level, supporting a lighthouse of weight F, 
the exact solution for the static vertical displacement is 

(22) u(x) = | (% 2 .x 2 (h(x) - - iPtax + F(& ~ X) + ab^ + fo^j Vx e [-a, b] 

where j t is the specific weight of the rock and fluid together, 7/ is the specific weight of the water, E is the 
modulus of elasticity for the rock, ipi and i/>2 are given as 

(23) ip 1 =j t + j j ip 2 = It - If 
and H(x) is the standard heaviside function 



(24) H{x) 



1 for x > 
for x < 



The exact solution was obtained by integrating the one dimensional boundary value problem for axial loading 
of a column, using the specific weights of water and the rock as body forces. Note that this solution is the 
exact solution of the set of differential equations that represent classical solid mechanics. It is not precisely 
the exact solution to the governing equations of the peridynamic formulation above; however, the results 
show good performance for the proposed formulation. 

The parameters used for this problem are given in Table [l] For completeness, the material constants 

have been presented in terms of Young's modulus, the bulk modulus, and the shear modulus. Any one of 
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Stress Pressure Displacement (mm) 

Figure 3. Lighthouse supported by a partially submerged column of rock: (left) problem 
geometry and applied stresses (right) verification results for the proposed formulation. The 
gravitational stress includes the weight of the lighthouse itself and the accumulated weight 
of the saturated rock column. 



Parameter 


Value 


a 


5 m 


b 


20 m 


F 


20 MN 


A 


78.54 m 2 


It 


25.67 kN/m 3 


If 


10.02 kN/m 3 


E 


29.7 GPa 


k 


9.0 GPa 




15.0 GPa 


7 


1.0 


6 


5.0 



Table 1. Lighthouse consolidation problem parameters 



these three can be calculated from the other two. The peridynamic analysis was performed for a 25 m long 
rectangular column of unit cross-sectional area discretized with elements of a volume 0.0015625 m 3 . The 
boundary conditions consist of fixing the bottom layer of elements such that the vertical displacement is 

10 



zero and applying a concentrated load to each element in the top layer consistent with the weight of the 
lighthouse distributed over the number of elements in the top layer. The gravitational force was applied via 
the body force term b(x) and the fluid pressure was applied via the pf term in equation (|9j). Figure [3] shows 
a comparison of the results for this problem using the proposed formulation and linear, one-dimensional, 
two-node bar finite elements. Although the error for both methods is relatively small, the finite element 
solution begins to deviate from the exact solution for the case that includes effective stress. 

The results illustrate that the modification of the pressure term according to equation ([9| does not 
invalidate the principle of correspondence. Even though in the proposed formulation the influence of the 
fluid pore pressure is engendered though the calculation of the force scalar state, the underlying material 
responds to loading in the same way as described by classical constitutive theory. 



4.2. Time dependent consolidation due to harmonic loading. For additional verification, we consider 
here the application of the proposed method to time dependent consolidation of a fluid saturated porous 
material under harmonic loading. This one dimensional, shown in Figure |4j consists of a column of material 
fixed at one end and loaded with an excitation force, F(t), at the other given by 

(25) F(t) = 50(1 - cos(75i)) 

The fluid pressure is assumed to vary according to the excitation force in time 



(26) p f (t) = 2F(t) ^1 

and linearly along the domain dropping from pt at the left side to zero at the right side where the excitation 
force is applied. These conditions correspond to a perfectly drained boundary at the surface of loading. 
At time r the pressure field reaches a state of zero pressure. The peridynamic analysis for this problem 
was performed on a unit-cross sectional domain of depth, L = 1 m discretized into cubic cells of volume 
0.0026653 m 3 . The problem parameters are listed in Table [2] The exact solution for the deflection of the 
surface on which the loading is applied can be obtained easily by integrating the one dimensional boundary 
value problem of classical elasticity and is given as 

(27) «*=o(*) = §(V(*) + ^ 

Figure [5] shows the harmonic loading and the maximum value of the fluid pressure, each measured at 
opposite ends of the domain. Figure [6] shows the evolution of surface deflection in time for both a dry 
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Figure 4. Problem geometry, boundary conditions and loading for the time dependent 
consolidation problem. 



Parameter 


Value 


F 


50(1 - cos(75i))kN 


T 


0.4 s 


E 


10.0 MPa 


k 


3.33 MPa 




5.0 MPa 


7 


1.0 


8 


3.5 



Table 2. Time dependent consolidation due to harmonic loading problem parameters 



porous material and one in which effective stress is taken into account. For both cases the proposed method 
performs well in capturing time dependent consolidation. 

4.3. Surface subsidence due to well extraction. The next problem demonstrates the surface deflection 
caused by a decrease in fluid pressure inside a reservoir as fluid is extracted. The drop in pore pressure from 
an initial state of equilibrium causes the porous material to contract volumetrically leading to subsidence. 
This problem is motivated by a similar problem studied in |13[ I18j . The fluid pressure field for this problem 
is based on the numerical results of these studies. A linear drop in pressure is assumed form the borehole 
radially to the outer boundary of the domain. The change in pressure from extraction is assumed to be -150 
psi at the borehole which represents an extraction duration of 75 days. The authors of [13l [18] report that 
the surface subsidence should be 4.5 ft after 75 days. 
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t (days) 

Figure 5. Evolution of the excitation force and fluid pressure in time for the time dependent 
consolidation due to harmonic loading problem. The fluid pressure is measured at the 
maximum depth of the domain (x = 10 m) and the excitation force is measured at the top 
surface of the domain (x = m) . 
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Figure 6. Evolution of surface deflection for the time dependent consolidation due to 
harmonic loading problem comparing the behavior including the effect of pore pressure and 
the behavior of a dry material. 

The peridynamic analysis of this problem using the proposed method was performed on a domain of 
dimensions shown in Figure [7] The field-scale do main was discretized with cells of inner radius of 40 ft 
and the parameters used are reported in Table [3j The boundary conditions consist of symmetry conditions 
for displacement along the outer edges of the domain and zero displacement in the X3 direction along the 



bottom surface. From the results shown in Figure [8] it is clear that the proposed formulation captures the 
same effective subsidence of the surface due to extraction as reported in previous studies. 



u X2 = 




Figure 7. Problem geometry and boundary conditions for the surface subsidence due to 
well extraction problem (all dimensions are in ft). 




Figure 8. Surface subsidence due to well extraction. Extraction from the reservoir leads 
to a decrease in fluid pressure near the borehole. This in turn results in consolidation of the 
reservoir material. The displacement in the a^-direction is measured along the o^-axis at 
the top surface of the domain. The loading conditions correspond to an extraction duration 
of 75 days. 



4.4. The effect of leak-off during hydraulic fracturing. With regard to applications in hydraulic 
fracturing, the present formulation is well suited for capturing the influence of leak-off effects. As fracturing 
fluid is injected into the reservoir, depending on the properties of the fluid and the injection rate, the fluid will 
penetrate into the geologic material to a varying extent. For high injection rates or geologic materials with 
low porosity very little leak-off occurs and the primary driver for damage propagation comes directly from 
the pressure of the fluid on the surface of the borehole. When the injection rate is lower or the penetrability 
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Parameter Value 

k 3.8642 ksi 

fj, 3.846.2 ksi 

Ap -150.0 psi 

7 1.0 

6 3.5 



Table 3. Surface subsidence due to well extraction problem parameters 



Parameter 


Value 


k 


9.0 GPa 




15.0 GPa 


7 


1.0 


S 


3.5 



Table 4. Hydraulic fracturing problem parameters 

of the material itself is high the fluid leaks into the pores of the geologic material leading to an increased state 
of stress in the porous media. To avoid the complexity of the interaction between crack growth and fluid 
penetration we only consider here the process for pristine material prior to crack initiation. The effective 
influence of leak-off on the porous material is measured by the dilation of borehole in terms of displacement 
for various assumed depths of fluid penetration. 

The geometry of the leak-off problem is shown in Figure [9] The domain was discretized with cells of 
volume 0.001153 cm 3 . The analysis was performed on a three dimensional domain with a thickness of 0.5 
cm assuming plane-strain conditions (no strain in or out of the plane of the paper) . Symmetry was used to 
reduce the problem to model only one quarter of the material in the vicinity of the borehole. The loading 
was prescribed as a constant force applied to each element equivalent to 1.0 Mpa of confining pressure on 
the top and right surface. On the borehole surface a force was prescribed on each element equivalent to 
a pressure of 0.1 Gpa. To preserve symmetry, no displacement was allowed in the ^-direction along the 
bottom surface and no displacement was allowed in the xi-direction along the left surface. The displacement 
was prescribed as zero in the a^-direction on the top and bottom surface. 

To investigate the effect of leak-off several simulations were run with varying degrees of fluid penetration. 
For each case, the fluid pressure was prescribed as 0.1 Gpa on the borehole surface and zero at the depth 
of penetration with a linear drop in pressure between. The results of this study are shown in Figure |10| 
which plots the displacement in the xi-direction along the xi-axis. To verify the solution comparisons are 
made, for the case of no fluid penetration with the results produced by Abaqus (a commercial finite element 
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5.0 cm 



1.0 cm 




4.0 cm 



Figure 9. Geometry of the leak-off problem. In the region of fluid penetration, the pressure 
is specified as a constant value at the borehole surface and zero at the outer boundary. 
Between these two surfaces, a linear drop in pressure is prescribed. There is no gradient in 
the pressure in the X3 direction. 

code) using eight-node brick elements with similar boundary conditions, loading, and element density. Good 
correspondence is obtained for this case between the two codes. 
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Figure 10. The effect of leak-off on the displacement in the vicinity of the borehole during 
hydraulic fracturing. The symbols show verification of the case without leak-off by compar- 
ison with Abaqus. The solid lines represent the effect of leak-off with varying depths (in 
centimeters) of fluid penetration beyond the borehole surface. 
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Borehole Edge Right Domain Boundary 
(xi = 1.0) (gi = 5.0) 



Dry 4.85 mm 

4.0 cm of leakoff 6.30 mm 



1.62 mm 
4.17 mm 



Table 5. Comparison of the displacements for the dry case and the case with 4.0 cm of fluid penetration 



Also shown in Figure 10 is the effect of varying the depth of fluid penetration. Near the borehole, 
the displacement is greatly increased due to leak-off effects for slight penetration, but the influence does 
not continue to grow for greater penetration depths. However, the response away from the borehole is 
almost unaffected by the initial penetration, but continues to increase rapidly as the fluid depth increases. 
Representative values of the displacement both at the borehole surface and the right edge of the domain are 
given in Table [5] for the dry case and the case of 4.0 cm of fluid penetration. 

If the pore fluid pressure is large enough, in comparison to the confinement pressure due to in-situ 
conditions, the presence of the fluid leads to increased stress in the porous material rather than stress relief (as 
was seen in the lighthouse consolidation example) and becomes the driving force behind damage propagation. 
Cases like this represent a fundamental deviation from the context of consolidation or subsidence problems 
for which the traditional effective stress principle was developed. 



5. SUMMARY 

This work presents a nonlocal model for fluid-structure interaction based on an effective stress concept. 
Through the principle of correspondence the underlying constitutive behavior of the porous material is 
unaltered, but the effect of fluid pore pressure is accurately captured. The method was demonstrated for a 
number of numerical examples including, one-dimensional consolidation and three-dimensional subsidence, 
the proposed formulation was also used to show the effect of leak-off during the pre-crack initiation stages 
of hydraulic fracture. The results of each of these examples show strong correlation with results previously 
reported and suggest that the formulation accurately captures the effective stress principle. This outcome is 
significant from the standpoint that it enables the development of more advanced methods for fluid-structure 
interaction that incorporate all the advantages of a non-local formulation such as discontinuity capturing and 
long range forces. Ultimately, it is intended that this work lay the foundation for the study of more complex 
interaction caused by flow induced damage of porous materials which can be used to model engineering 
problems like hydraulic fracturing. 
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